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VLBI for Gravity Probe B: The Guide Star, IM Pegasi 

N. Barteb, M. F. Bietenholz^’^, D. E. Lebach^, R. R. Ransom^’^, M. I. Ratner^, and I. I. Shapiro^ 

ABSTRACT 

We review the radio very long baseline interferometry (VLBI) observations of the 
guide star, IM Peg, and three compact extragalactic reference sources, made in support 
of the NASA/Stanford gyroscope relativity mission, GP-B. The main goal of the ob¬ 
servations was the determination of the proper motion of IM Peg relative to the distant 
universe. VLBI observations made between 1997 and 2005 yield a proper motion of 
IM Peg of —20.83 ± 0.09 mas yr“^ in a and —27.27 ± 0.09 mas yr“^ in d in a celestial 
reference frame of extragalactic radio galaxies and quasars virtually identical to the In¬ 
ternational Celestial Reference Frame 2 (ICRF2). They also yield a parallax for IM Peg 
of 10.37 ± 0.07 mas, corresponding to a distance of 96.4 ± 0.7 pc. The uncertainties 
are standard errors with statistical and estimated systematic contributions added in 
quadrature. These results met the pre-launch requirements of the GP-B mission to not 
discernibly degrade the estimates of the geodetic and frame-dragging effects. 

Subject headings: astrometry — binaries: close — gravitation — radio continuum: 
stars — radio continuum: galaxies — relativity — stars: activity — stars: individual 
(IM Pegasi) — techniques: interferometric 


1. Introduction 

Einstein’s theory of general relativity (GR) predicts that space-time is curved and warped 
by the mass and the angular momentum of a gravitating body. The kinematics of objects in the 
vicinity of such a body are expected to differ from those based on Newton’s theory of gravity. 
The NASA/Stanford spaceborne mission Gravity Probe B (GP-B) was designed to measure these 
differences with four gyroscopes essentially freely falling around Earth, in a polar orbit at an altitude 
of about 640 km. According to GR, the gyroscopes’ spin axes should process by ~6"6 yr ^ in a 
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north-south direction due to the curving of space-time by the Earth’s mass, and by ~39 mas yr“^ 
in an eastward direction due to the warping of space-time by the Earth’s angular momentum. 
These precessions are dubbed the “geodetic” effect and the “frame-dragging,” (“Lense-Thirring,” 
or “gravitomagnetic”) effect, respectively, and can be considered as rotations of any near-Earth 
inertial frame relative to the distant universe. 


The gyroscopes were set spinning within a quartz block with respect to which, through su¬ 
perconducting quantum interference devices (SQUIDs), each of the spin axes’ precessions could 
be measured with an accuracy of 100 fias within ~40 days. The quartz block was bonded to a 
~15 cm diameter quartz optical telescope. All of these devices — the quartz block, the SQUIDs, 
and the quartz telescope (together termed the probe) — were placed in a dewar filled with liquid 
helium and cooled to a temperature of 1.8 K. Ideally, the orientation of the probe could be best 
provided by pointing the telescope to one suitable quasar, one of the most distant compact types 
of single object in the universe. However, quasars were too dim for the telescope to track. Only 
stars in our galaxy were bright enough. The challenge was therefore to hnd a “guide star,” bright 
enough for the on-board telescope to track, sufficiently isolated to limit light contamination, and 
suitably located to not unnecessarily decrease the accuracy of the measurement of the precessions 
of the gyroscopes. Then the guide star’s motion on the sky needed to be measured with respect 
to quasars or distant galaxies so that the precession of the gyroscopes could be determined with 
respect to the distant universe. 


We described the astronomical effort undertaken with the guide sta r for the support of the 
GP-B mission, “VLBI for Gravity Probe B,” in seven papers: “Overview” (jShapiro et al.ll2012l . Pa- 
per I); “Monitoring of the structure of the reference sources 3C 454.3, B2250-I-194, and B2252-I-172” 
( Ransom et al.ll2012al . Paper II); “A limit on the proper motion of the ‘core’ of the quasar 3C 454.3” 


suits from other technio 

ues” ( 

iiebach et al. 2012. Paoer IVj; “Proper motion and parallax of the 

guide star, IM Pegasi” ( 

Ratner et al. 20121. Paoer Vj: “The orbit of IM Peeasi and the location of 

the source of radio emission” 

(Ransom et a.l. 

2012b. Paper VD; and “The evolution of the radio 

structure of IM Pegasi” 

Bietenholz et al. 2012 

, Paper VII). Here we give a synopsis of these papers. 


2. Search for the best guide star and selection of IM Peg 

A suitable guide star for GP-B had to simultaneously meet a number of requirements: 

1. The motion of the star needed to be known or measurable with an uncertainty small enough 
not to significantly increase the prelaunch anticipated standard error of GP-B of <0.5 mas yr“^ 
in each of the two precessions. The specific requirement was that the standard error be 
<0.14 mas yr“^ for measurements of the star’s motion in each coordinate. In 1990, when we 
entered the hrst stage of our search for a suitable guide star, no star’s (or any other celestial 
object’s) motion on the sky was known with such accuracy in the optical. Such accuracy 


































- 3 - 


had previously been achieved only in the radio. Pioneered bv IShaoiro et ahl (jl979|), astro- 
metric measurements with the radio technique of very long baseline interferometry (VLBI) 
of a quasar relative to another quasar nearby on the sky yielded a proper-motion estimate 
with a standard e rror less than half of the requirement for GP-B (jMarcaide &: ShaDirolll983l : 


Bartel et al.lll986l b Therefore VLBI appeared to be at that time the only option to determine 


the motion of a suitable guide star with the desired accuracy. 

2. The star needed to be sufficiently bright in the radio, preferably with a flux density at least 
of order 1 mJy, for VLBI astrometry measurements to be successful. 

3. The star needed to be sufficiently bright in the optical, at least 7th magnitude, for the on¬ 
board telescope to detect. 


4. The star’s angular distance, D, from the ecliptic needed to be between 20° and 40°. The lower 
limit was set by the requirement that, with a suitable Sun shield at the telescope, sunlight 
not enter the dewar and boil off the liquid helium. The upper limit was set by the power 
requirements of the spacecraft, since the solar panels were fixed to the spacecraft, which would 
continually point to the guide star, making the power received by the solar panels dependent 
on D. 


5. The star needed to be located as close as possible to the celestial equator, so as to maximize the 
sensitivity of the measurement of the frame-dragging precession, which essentially decreases 
with the cosine of the star’s declination. 

6. The star needed to be sufficiently isolated on the sky. Neighbouring stars or reflection nebulae 
could possibly intolerably decrease the accuracy of the pointing of the telescope to the star’s 
center of brightness. 


To find the most suitable guide star for GP-B, we conducted a VLA survey at 8.4 GHz from 
1990 to 1992 of ~1200 stars with V magnitude <6.0 and a declination between —20° and -|-40°. No 
single stars were found in the survey with a flux density above ~1 mJy. In fact, we detected only 
previously k nown radio stars, all of them binaries of the RS Canum Venaticorum type (see also, for a 
later survey. iHelfand et al.lll999l ). We identified four potentially suitable guide stars: A Andromeda 
(HR 8961, -|-46° declination), HR 1099 (-1-1°), HR 5110 (-1-37°), and IM Peg (HR 8703, -|-17°). 
After investigating each of the four candidates with VLBI, studying the sky fields around them in 
the optical, and checking on possible VLBI reference sources for each of them, we selected IM Peg 
as the most suitable for GP-B. 


2.1. Properties of IM Peg and its surroundings 

IM Peg is a binary with a giant primary and Sun-type secondary. We summarize its optical 
characteristics and previously known properties in Tabled! The radius of the primary is ~13.3 ± 
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0.6 Rq, which translates at a distance of IM Peg of ~100 pc, to an angular radius of ~0.64 ± 0.03 
mas. This radius is comparable to the full-width at half-maximum (FWHM) of the synthesized 
beam of a global VLBI array operating at 8.4 GHz, and therefore in principle allows for an accurate 
determination of the location of the radio emission with respect to the center of the primary. 


A prerequisite for successful guide star tracking was that the onboard telescope be able to lock 
sufficiently accurately on the center of the optical disk of the star (or a point with an offset from this 
center constant during the mission). One concern here was the appearance of variable dark spots 
on the primary, shown by Doppler imaging and photometry to cover ~15% of the star’s surface. 
If the variations in the size or location of these spots resulted in a linear trend over the course of 
the GP-B mission in the centroid of optical brightness, relative to the center of the disk of the IM 
Peg primary, then that trend would have the same effect on the GP-B data reduction as an equal 
error in our VLBI measurement of the proper motion of IM Peg. Such a trend is not necessarily 
implausible, given the year-to-year variations in IM Peg’s optical brightness and the astrophysical 
plausibility of variations in typic al spot latitudes analog ous to the systematic variations of sunspot 


locations on the Sun. However, iMarsden et al.l (120071 ) showed, through extensive spot mapping 


spanning the duration of the GP-B mission, that any error due to variability in the spot pattern 
would yield no more than a 0.04 mas yr“^ drift in the optical center. 


A second concern was the surroundings of IM Peg on the sky. Any star or nebula within the 
field of view of the telescope when locked on IM Peg, possibly combined with IM Peg’s photometric 
variability, could cause systematic astrometric errors. However, for IM Peg, no other star within, 
e.g. 12', is brighter than V magnitude 10, which can be compared with the corresponding relative 
large brightness of IM Peg, that varied during the mission only between V magnitude 5.7 and 6.0. 
All in all, through extensive optical and even millimetric observations of the CO(J=l—0) line (the 
latter to search for any molecular cloud that could be associated with a reflection nebula), we 
constrained any astrometric systematic error to a negligible value: The combined error from dark 
spots, neighbouring stars, and a nebula, combined with the photometric variability of IM Peg, we 
estimated to be smaller than 0.05 mas yr“^. 


3. Radio reference sources for IM Peg in the distant universe 

Our goal was to determine the motion on the sky of IM Peg relative to the distant universe, 
with a standard error not exceeding 0.14 mas yr“^ in either coordinate. To achieve this goal, we 
selected radio reference sources with respect to which the motion of IM Peg could be measured. 
These sources needed to be at cosmological distances, be sufficiently strong and compact for VLBI 
observations, and located nearby to IM Peg on the sky so as to minimize systematic astrometric 
errors. 

We selected three reference sources, with the strong quasar, 3C454.3, as the main one. The 
other two were the active nucleus of the radio galaxy, B2250-I-194, and the unidentified source 
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Table 1. Optical characteristics and previously known properties of IM Peg 


Parameter 

Value 

Reference 

Hipparcos, q;(J 2000), epoch 1991.25 

22*1 53™ 02?278706 ± 0.63 mas 

1 

Hipparcos, (5(J2000), epoch 1991.25 

16° 50' 28'.'53982 ± 0.43 mas 

1 

Hipparcos, parallax (mas) 

10.33 ±0.76 

1 

Hipparcos, distance (pc) 

96.8 

+ 7.7 
- 6.2 


Stellar Properties^ 

V magnitude range during mission 

5.7 to 6.0 


2 

Mass (Mq) 

1.8 ±0.2 

1.0±0.1 

3,3 

Spectral Type 

K2 III 

G V*' 

4 

Teff (K) 

4550 ± 50 

5650 ± 200^' 

4,3 

Radius (R©) 

13.3 ±0.6 

1.00 ±0.07*' 

4,3 

Radius (mas)'^ 

0.64 ± 0.03 

0.048 ± 0.004‘> 

4,3 

Orbital Elements'^ 

a sin 2 (R©) 

16.70 ±0.02 

30.34 ±0.03 

3,3 

a sin 2 (mas)'^ 

0.806 

1.464 


P (days) 

24.64877 ± 0.00003 

3 

i n 

65 to 80, >55 

4,5 

e 

0.0 (assumed) 

4 

Tconj (HJD)d 

2450342.905 ±0.004 

3 


^Two entries correspond to the two stars of the binary system, with entries for the primary 
listed first. First reference is for the first entry, second reference, if present, is for the second 
entry. 

^The spectral type, effective temperature, and radius of the secondary are inferred from the 
flux ratios (at two wavelengths) of the two stellar components and the values for the radius 
and effective temperature of the primary under the assumption that the secondary is a main 
sequence star. 

‘^Computed for a system distance of 96.4 di 0.7 pc. The uncertainty in the asini value in Rq 
units is not propagated into mas, since the uncertainty in the inclination is the dominant source 
of error in any spectroscopic determination of the semimajor axis. 

'^Heliocentric time of conjunction with the K2 III primary behind the secondary. 


References. — 1 . Hip parcos Catalogue teSAlll997|h 2. Pape r V: 3. iMarsden et al] ll2005l~l : 4. 
iBerdvugina et ahl lll99^ j (e = 0.006 ib 0.007); 5. Ibebach et al.l lll999|j . 
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(most likely also a quasar or a radio galaxy), B2252+172. We display their positions on the sky, 
together with that of IM Peg, in Figure [H and list their characteristics in Table [2j In Figure [2] we 
present typical VLBI images at 8.4 GHz of IM Peg and each of the three references sources. 


All the sources are approximately located along the same axis (north-south), allowing perhaps 
for a more accurate relative correction of tropospheric and ionospheric effects. However, each of 
the three sources has advantages and disadvantages as reference sources for IM Peg. The quasar 
3C 454.3 is both the strongest a nd the closest on the sky to IM Peg. It was used as a reference source 
for IM Peg from 1991 to 1994 ([Lestrade et al.lll999l b and could therefore extend the time baseline 
of our VLBI observations. The disadvantage is the complexity of 3C 454.3’s structure and its 
superluminal nature. In Figure[2]we identified six components of which the easternmost component. 
Cl, is still compact at the highest angular resolution, has a flat or inverted high-frequency spectrum 
(jPagels et al.l 120041 ). and is located almost exactly a t the same position (Paper H) as the core of 
the quasar seen at 43 GHz (|Jorstad et al.ll200ll . l2005l L The component Cl was therefore considered 
to be the core at 8.4 GHz, the part most closely related to the supermassive black hole and to 
the center of mass of the quasar; it was thus the reference point in the quasar’s structure for our 
astrometric measurements. The other components are moving away from Cl, with apparent speeds 
of up to ~5 c (Paper H). The remaining sources were both rather compact, so that the brightness 
peak could simply be taken as the reference point in the structure; their disadvantages were that 
B2250-I-194 was located relatively far away from IM Peg and B2252-I-172 was rather faint (with an 
unknown redshift). 


In addition, we used geodetic and astrometric VLBI observations regularly conducted over 
many years to determine the limit on the proper motion of two of the three reference sources, 
3C 454.3 and B2250-I-194. These determinations were made without consideration of a particular 
reference point in the structure of the two sources but were made with respect to a celestial reference 
frame, CRF. This frame is defined by a large number of extragalactic radio sources distributed over 
the sky and is for our purposes essentially the same as the International Celestial Reference Frame, 


Table 2. Characteristics of the sources 


Source 

Type 

Separation 
Ao(°) A6{°) 

Flux density® 

(Jy) 

Redshift 

Distance ^ 
(Mpc) 

3C 454.3 

quasar 



7-10 

0.859 

1610 

B2250-I-194 

galaxy 

-0.2 

3.6 

0.35 - 0.45 

0.28 

880 

B2252-I-172 

unidentified 

0.4 

1.4 

0.017 



IM Peg 

RS CVn 

-0.1 

0.7 

0.0002 - 0.08 

0.0 

0.0 


^The range gives the lowest and highest flux density measured at 8.4 GHz with the VLA 
during the course of the observations, 1997 January to 2005 July. 

’^The angular diameter distance for a flat universe with Hubble constant, 
Hq=7{) km s~^ Mpc“^, and normalized density parameters, Qm = 0.27 and $2^=0.73 . 
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Fig. 1.— Relative positions (J2000) on the sky of the four radio sources used for GP-B astrometry. 
The a (east-west) and 5 (north-south) directions on the plot are on the same scale. Figure taken 
from Paper III. 




MilMARC SEC 
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Fig. 2.— VLBI images at 8.4 GHz of IM Peg (upper left panel) on 2004 December 11, and the three 
reference sources, 3C 454.3 (lower left panel) also on 2004 December 11, B2250+194 (upper right 
panel) on 2000 November 5, and B2252+172 (lower right panel) on 2003 May 18. The restoring 
beam characterizing the resolution in the image is given for each panel at its 50% contour (FWHM) 
in the inset. North is up and east to the left. For more information, see Paper II and Paper VII. 
Figure taken from Paper I, but rearranged. 




- 9 - 


ICRF2 ( Fev et al.l 120091) 


All in all, we took Cl in 3C 454.3 as the primary reference point with respect to which the 
motion of IM Peg was measured. The other measurements served to determine the stability of Cl 
on the sky so that finally the motion of IM Peg could be determined with respect to the distant 
universe. 


4. Astrometric VLBI observations of IM Peg and its radio reference sonrces, and 

data analysis 

To determine the motion of IM Peg relative to Cl and test Cl’s stability relative to the 
brightness peaks of the other two reference sources, we took into account: 

1. IM Peg is catalogued as a binary. However, we could not rule out that it is perhaps part 
of an extended triple or even larger multiple system. Therefore, to ensure our obtaining a 
sufficiently accurate value of the proper motion of IM Peg for the ~15 months GP-B data- 
collection interval in 2004 and 2005, we spread a sufficiently large number of VLBI measure¬ 
ments over eight and a half years to be able to determine any nearly constant change of the 
proper motion with time (“proper acceleration”), or place a limit on it. 

2. For the best determination of IM Peg’s parallax, we spread the VLBI measurement epochs 
appropriately to cover the different seasons each year. Similarly, to allow us to model any 
orbital component of the motion of the radio emission, we carefully distributed the epochs 
over the phase of the binary orbit, which could be computed from earlier spectroscopic results. 

3. To measure the time-variable brightness distributions of IM Peg and the three reference 
sources, and in particular to identify Cl in 3C 454.3, we selected each VLBI observation 
session to be long enough to ensure that the u-v coverage was sufficiently dense for us to 
obtain excellent imaging of the sources. 

4. Since we expected IM Peg to fluctuate strongly in radio brightness, and since we knew that 
radio emission from B2252-I-174 was relatively weak, we needed our VLBI array to be suffi¬ 
ciently sensitive to nearly guarantee detection of the sources in each observing session. 

5. Lastly, we needed an observing frequency for which receivers were available at each of the 
telescopes of the VLBI array, and which guaranteed sufficiently high angular resolution and 
low system temperatures. 

All said, 35 sessions of astrometric VLBI observations were made between 1997 January 16 and 
2005 July 16 with about four sessions every year, each about 11 to 15 h long. The observations were 
made with a global array of 12 to 16 radio telescopes comprised of most or all of: MPIfR’s 100 m 
telescope at Effelsberg, Germany; NASA/Caltech/JPL’s 70 m DSN telescopes at Robledo, Spain, 
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Goldstone, CA, and Tidbinbilla, Australia; NRAO’s ten 25 m telescopes of the VLBA, across the 
U.S.A; NRAO’s phased VLA, equivalent to a 130 m telescope, near Socorro, NM; and, at early 
times, NRCan’s 46 m Algonquin Radio Telescope near Pembroke, ONT, Canada, and NRAO’s 
43 m telescope in Green Bank, WV. The observing frequency was 8.4 GHz. 


All of these VLBI observations were made to extract interferometric phase delays related to 
the difference of arrival times of a radio wave from a celestial source at each pair of antennas of a 


VLB I array. Such observations allow the most accurate astrometric measurements (|ShaDiro et al 


1979|) and can yield relative positions and proper motions referenced to particular points in the 
structure of the sources with uncertainties as low as ~10 /ras and ~10 /ras yr“^, respectively (for 
the earliest such measurements, see, e.g., Marcaide &: Shapiro 1983; Bartel et al. 1986; and for a 
recent review, see Reid &: Honma 2014). 


The phase delay is given by 

( 1 ) 

where ((>(w, t) is the interferometric phase or “phase” at observing angular frequency, uj, and time, 
t. The phase can be expressed as 


(j)(^iO,t') — Ul [Tgeom(0 T Rnst(^)^) T T ’^struc(^)^) T '7'noise(^) ^)] T 27rA^(w,t), (2) 

where Tgeom(i) is the “geometric delay,” the difference in the arrival times of the radio wave in 
vacuum at the two antennas; Tinst(w, t) describes the difference in the instrumental delays (including 
clock behavior) at the two antenna sites; ratm(w,t) gives the difference in radio wave propagation 
times to the two antennas due to all tropospheric and ionospheric effects; Tstruc{^,t) is the delay 
contribution from source structure to account for any non-pointlike brightness distribution of the 
celestial source; rnoise(‘^j0 is the (thermal) noise contribution to the phase measurement; and 
N{u}, t) describes the integer number of 27r ambiguities, or “phase wraps,” in the measurement. 

The astrometric information of a VLBI measurement is given by the geometric delay which, 
apart from relativistic contributions and those due to Earth’s motion relative to the solar system 
barycenter, is given by 

rgeom(i) = - [Bit) ■ s(t)] , (3) 

c 

where c is the speed of light in vacuum, B(t) is the 3-dimensional vector between two antennas of 
a baseline of a VLBI array, and s(t) is the unit vector in the direction of the observed source. 

For phase-delay VLBI observations to succeed, to be able to remove the 27r phase ambiguities, 
and to reduce many sources of astrometric errors, we needed to switch the antennas rapidly between 
sources. For the first 23 sessions, we used a typical sequence of 3G 454.3 (80 s) - IM Peg (170 s) - 
B2250-I-194 (80 s). For the remaining 12 sessions, B2252-I-172 was included in the observations and 
the sequence was altered to 3C 454.3 (80 s) - IM Peg (125 s) - B2250-I-194 (80 s) - 3C 454.3 (80 s) 
- IM Peg (125 s) - B2250-I-194 (80 s) - B2252-I-172 (90 s). In addition to the astrometric 8.4 GHz 
observations, we also observed once at 5.0 GHz and once at 15.0 GHz, to allow investigation of 
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the compactness and spectral properties of the components of the sources. Almost all observations 
were recorded in both right and left circular polarization, and all were processed with the VLBA 
correlator at Socorro, NM. 

We used a custom software package (Paper IV) that allowed us to analyze the data efficiently. 
First, we connected the phases and removed the 2ti phase ambiguities to generate the correct phase 
delays. Then we corrected the phase delays for the effects of the troposphere and the ionosphere 
with a priori estimates. We further processed them with a Kalman hlter to model the residuals 
of the troposphere and the ionosphere, and the clock offset at each telescope from the one chosen 
as the reference clock for the VLBI array. Then, for the more extended source 3C 454.3, the 
structure effects were removed from the phase delays based on the reference point. Cl, defined 
with images made with NRAO’s AIPS software package. The other two reference sources were 
sufficiently compact that no correction for their source structure was needed. 

The relative weakness of IM Peg’s radio emission, and sometimes also that of B2252+172, 
made it necessary for us to develop a special technique for the analysis of GP-B VLBI data. It 
combines the advantages of parametric model-fitting via weighted least-squares applicable for the 
relatively strong sources 3C 454.3 and B2250-I-194 with the sensitivity of phase-referenced mapping, 
the latter required for IM Peg and B2252-I-172. The merged analysis technique yielded superior 
results to the phase-referenced mapping technique (Paper IV). 

The radio emission from IM Peg varied from session to session, reflecting the complex astro- 
physical nature of the environment of the star. In most cases the images of IM Peg showed only 
one clearly defined component. In these cases the coordinates were determined for the maximum 
of a two-dimensional Gaussian fit to that component. In the remaining cases, there was more than 
one local maximum of approximately equal brightness in the emission region: in eight cases there 
were two such maxima, and in one case there were three, all of them significant. In such cases, we 
took the position of IM Peg to be the mean position of the different maxima. 

The final result of the analysis of the GP-B VLBI observations was a set of coordinates with 
statistical standard errors for each source and each observing epoch. These sets were the basis for 
the study of the motion, or upper bounds on it, of the sources with respect to each other. 

In addition, we analyzed geodetic astrometric VLBI observations and made use of interfero¬ 
metric group delays. Group delays are given by 

doj 

These are observations of thousands of sources from all over the sky. These observations are done 
on a routine basis, independent of the GP-B program. For 3G 454.3 we used all available data 
from a total of 1,119 observing sessions, from 1980 to 2008. In support of GP-B, the second 
reference source, B2250-I-194, was included in 38 sessions of routine geodetic astrometric group- 
delay observations between 1997 and 2008. Group-delay observations have the advantage that 
27r ambiguities are not inherent in the data and that the positions are determined relative to a 
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celestial reference frame (CRF). The disadvantage, however, is that the parameter, rstruc(‘^) i); was 
not determined and that therefore the exact reference point in the brightness distribution of the 
sources for the position measurement remained undefined. 


All available geodetic astrometric data with a total of 6.5 million group-delay determinations 
were processed via least-squares and yielded estimates of coordinates with statistical standard errors 
for 3C 454.3 and B2250-I-194 for each observing session, keeping the coordinates of the other sources 
constant. This set of solutio ns forms our C RF, which for G P-B^s purposes is virtually identical 
(see, e.g., iPetrov et al.l 1200911 to the ICRF2 (jFev et al.N2009l b the most fundamental such frame 
presently in use. 


5. Positional stability of IM Peg’s reference sources 

The phase-delay determinations of the position of Cl, the core of 3C 454.3, relative to those 
of the brightness peaks of B2250-I-194 and B2252-I-172, are plotted in Figure El The group-delay 
determinations of the position of B2250-I-194 relative to the CRF are also plotted in Figure El See 
Paper III for a similar plot of the group-delay determinations of the position of 3C 454.3 relative 
to the CRF. 

Weighted least-squares straight line fits to each set of position determinations yielded the 
position at epoch and proper motion for each coordinate separately, along with their statistical 
standard errors and correlations. In each case the fits to the data gave Xu ix^ degree of 
freedom) larger than unity, indicating systematic errors in the position determinations. These 
errors can for instance be due to uncorrected effects from source structure and/or atmospheric 
and ionospheric variations, and are difficult to quantify. Therefore the total standard errors of the 
position determinations were determined partly empirically by including a particular constant for 
each fit. This constant was added in quadrature to the statistical standard errors for each set of 
position determinations and for each coordinate separately so as to obtain = 1 for each fit. 

The resultant proper-motion values with their standard errors are listed for each set of position 
determinations in the panels of Figure El The upper panel of Figure El indicates that there appears 
to be some motion of Cl relative to B2250-I-194 at the 3.5a significance level. The source 3C 454.3 
has a long jet, indicating strong activity. Some motion of a component like Cl even if located close 
to the black hole could therefore be expected. No significant motion is discernible for B2250-I-194 
relative to B2252-I-172. No matter whether there is some motion due to activity near the black 
hole or not, the la upper limits of the proper-motion values for Cl relative to B2250-I-194 and 
B2250-I-194 relative to the CRF are each smaller than 40 //as yr“^. Combining the phase-delay 
proper-motion determination of Cl relative to B2250-I-194 with the group-delay proper-motion 
determination of B2250-I-194 relative to the CRF gives for the time period from 1998 to 2005 a 
Icr standard error for the mean position of Cl in the CRF of 45 and 68 //as. The proper motion 
of Cl relative to the CRF for that same period is 33 ± 13 //as yr“^ and —35 ± 21 //as yr“^ in a 
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Fig. 3.— Top panel: The estimated coordinates of Cl, the core in 3C 454.3, relative to those 
of the brightness peaks in B2250+194 and B2252+172 (except for offsets) as a function of time 
from GP-B VLBI phase-delay observations. The lines give the fit to these coordinates. The values 
for the slopes of the lines and the statistical standard errors for = 1 are also given. Bottom 
panel: The relative coordinates of B2250-I-194 as determined from routine geodetic group-delay 
VLBI observations of up to ~4000 extragalactic sources scattered over the sky. Figure taken from 
Paper III. 





-14- 


and S, respectively, which we do not consider significant. We therefore give a la upper limit on 
the magnitude of the proper-motion components of 46 and 56 / 2 as yr~^, respectively. The latter 
pair of determinations is our limit on the level of stationarity of the reference point Cl in 3C 454.3 
relative to the distant universe (Paper III). 


6. Motions of IM Peg relative to the distant universe 

With the stationarity of the reference point in 3C 454.3 determined, the motion of IM Peg 
relative to the distant universe could be analyzed. Five kinds of motion needed to be considered: 

1. proper motion mainly due to the motion of IM Peg relative to the motion of the Sun in the 
Galaxy, 

2. annual parallax due to the orbital motion of Earth around the Sun, together with the much 
smaller contribution of the motion of the Sun with respect to the solar system barycenter, 

3. change of proper motion over time due to a hypothetical third member of the IM Peg binary 
system, 

4. orbital motion around the common center of mass of the binary, and 

5. erratic motion of the center of radio emission within the binary system. 


6.1. Astrometric solutions 


To determine the proper motion, parallax, proper acceleration, and orbital motion, we used 
weighted least-squares to fit a linearized model to the 35 determinations of the positions of the 
reference point in IM Peg relative to the position of Cl in 3C 454.3. The uncertainties of these po¬ 
sition determinations are hard to compute theoretically. We therefore took them from an empirical 
analysis based on the computation of the root-mean-square (rms) scatter separately in a and 6, of 
the postfit residuals obtained for the 12 positions of B2252-I-172 relative to Cl and the 35 positions 
of B2250-I-194 also relative to Cl. The resulting uncertainty was 0.06 mas in each coordinate. In 
addition there is the previously discussed uncertainty of the position of Cl and the upper limit of 
the proper motion of Cl in the CRF that we considered for the standard errors of the position 
determinations of IM Peg relative to the CRF. 


The parameters of the model fit were the position of IM Peg at epoch, proper motion, parallax, 
and four sca lar parameters for the d escription of the projection on the sky of an orbit with zero- 
eccentricitv (IBerdvugina et al.lll999l ) and a period known from optical spectroscopic observations 


(Marsden et al. 


rdvupi 

l2nn.5l l. 


The fit resulted in rms values of the residuals somewhat different for each 
coordinate but much larger than the standard errors of the position determinations of IM Peg. The 
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source of this large scatter is the astr ophysical nature of t he radio emission, characterized by motion 
related to radio brightness changes ([Lebach et al.lll999ll and plausib ly in part to stellm ma gnetic 
field changes implied by spot maps from optical spectroscopy le.g.. iBerdyugina et alJl2000l '). For 
the fit we used uniform weighting independently for a and 5, and allowed their errors to haye non¬ 
zero correlation. After iteration, we obtained the final set of parameter estimates with statistical 
standard errors and = 1. This set is presented in Table El 

For the determination of the proper acceleration, we enlarged the set of positions of 35 epochs 
by four additional positio ns of IM Peg obtained at epochs between 1991 and 1994, also at 8.4 GHz, 
by iLestrade et al.l (|l995l l in support of the Hipparcos mission. These observations considerably 
extended our time baseline; howeyer, they were only of limited use for the estimate of parameters 
other than the proper acceleration because of inferior u-v coverage, lower angular resolution, and 
the position determinations not being referred to Cl. We found no significant proper acceleration; 
i.e., none larger than one statistical standard error; and so we include no such proper acceleration 
in the fit used to obtain our final results in Table El 


Figure m displays the 39 position determinations and the fit to the 35 position determinations, 
with the four earliest positions not used in the fit. The effect of the proper motion and the parallax 
can be clearly seen in the data as well as in the fit. Figure El displays the effect of the parallax 
still more clearly. It shows the 35 position determinations with the estimated position at epoch, 
proper motion, and the primary’s orbital motion subtracted. It also shows the estimated parallax 
ellipse. Figure El displays the orbit clearly. Again we plot the 35 position determinations with 
the estimated position at epoch, proper motion, but now with parallax rather than orbital motion 
subtracted. Also shown are the estimated orbit and the corresponding positions predicted from the 
orbit model. Figure [71 shows the 39 position determinations and the fit to the above 35 positions, 
but now with all model contributions, namely the estimated position at epoch, proper motion, 
parallax, and the primary’s orbital motion subtracted. We compare them to the corresponding flux 
densities of IM Peg. The scatter of the position residuals appears to be random with certainly no 
correlation to the measured flux densities. The rms of the scatter is ~0.4 mas in each coordinate, 
much larger than the standard errors of the position determinations and most likely dominated by 
the fluctuation of the positions of the stellar radio emission relative to the center of the primary 
of IM Peg. Figure El gives an artist’s three-dimensional rendition of IM Peg with the primary as a 
giant with dark spots and the secondary as a sun-like star, each in its estimated orbit. 


6.2. Systematic and total errors 

A detailed analysis of systematic errors affecting each of the nine fit parameters of IM Peg’s 
motion relative to the distant universe led us to the following conclusions. For the position at 
the mean epoch of the VLBI observations of 2001.29, the systematic error in each coordinate was 
computed from 0.5 times the angular radius of the primary along the sky projection of the normal 
to the binary orbital plane (p.a. = 130.5° ± 8.6°). For the later epoch of 2005.08 in Table El the 
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Fig. 4.— The 35 position estimates of IM Peg relative to Cl. The coordinate origin is arbitrary. 
The 2a errors (here and hereafter in the figure captions in the sense of the full length of ±la error 
bars) are much smaller than the symbols. The curve shows the fit of the solution, with orbital 
terms omitted for clarity. In addition, the four positions from 1991 to 1994 are shown (upper 
left); however, they were not used in the fit. The arrow indicates the direction and its length the 
magnitude of the proper motion in one year. Figure taken from Paper V. 
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Aa cos((5) (mas) 

Fig. 5.— The parallax ellipse (solid line) derived from the nine-parameter fit to the set of 35 
position determinations of IM Peg. The open circles give the observed position for each epoch after 
subtraction of the estimated position at the reference epoch, proper motion, and orbital motion. 
The 2 (t errors (of position determinations relative to Cl) are about a third of the size of the circles. 
The relatively large scatter is most likely dominated by the fluctuations of the positions of the 
stellar radio emission relative to the center of the disk of the IM Peg primary. The model positions 
(not shown) are almost exactly on the plotted parallax ellipse and close to the observed positions. 
The parallax ellipse, however, is slightly time-dependent due to the influence of planets. Over the 
course of our observations of several years, the parallax ellipse in the figure shifts by about the 
width of the line. 
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Fig. 6.— The orbit (solid line) of the primary of IM Peg derived from the nine-parameter fit to 
the set of 35 position determinations of IM Peg with the arrow giving the direction of motion. The 
large triangle on the upper left segment of the ellipse marks the ascending node of the inferred 
orbit, under the convention that at an ascending node an orbiting stellar component recedes from 
Earth. Thus the part of the ellipse with the arrow on it is closest to Earth. The open circles give 
the observed position for each epoch after subtraction of the estimated position at the reference 
epoch, proper motion, and parallax. The 2a errors (for the position determinations relative to 
Cl) are about twice the size of the circles. A solid line connects each observed position with the 
corresponding fitted position indicated by a dot on the estimated orbit. The large scatter of the 
observed positions about the model orbit ellipse is again most likely dominated by the fluctuations 
of the positions of the stellar radio emission relative to the center of the disk of the IM Peg primary. 
Figure adapted from Paper VI. 
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1995 2000 2005 


Observation Epoch 

Fig. 7.— The upper panel shows flux densities for our 39 VLBI sessions, either as a range of values 
measured with the VLA within a session or as a single value with la error bars obtained mostly 
from the VLBI images. The two lower panels show the position residuals for all sessions. We plot 
unweighted points as open circles and weighted points as either closed circles or triangles, with each 
triangle indicating a position computed as the mean position of the two or three resolved peaks in 
the stellar radio image for that session. The 2cj errors (for position determinations relative to Cl) 
are about the size of the symbols. For more information, see Paper V, from which this figure was 
taken. 
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Fig. 8.— Artist’s three-dimensional rendering of the IM Peg binary system as seen from Earth. The 
primary is the larger red star with dark spots, with the spin axis parallel to the orbit normal, and 
mass ejection along magnetic field lines in the polar region that could give rise to the observed radio 
emission. The blue contours show, as an example, the radio brightness distribution as observed 
on 2003 May 18 relative to the primary (see Paper VII). Similar radio emission relative to the 
primary at approximately equal phases in its orbit were observed on 1997 November 30 and 2000 
August 7 as shown in Figure [H The secondary is the smaller yellow star. The projected orbit of the 
primary is the same as the inferred radio source orbit shown in Figured The projected orbit of the 
secondary and the diameters of the two stars correspond to the nominal values given in Table [TJ 
The system is shown with the primary at its ascending node. Note that the primary’s visible pole 
is near its southeast side. The size of the figure is approximately 3.8 mas by 3.4 mas. 
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value for each component was added in quadrature with a possible maximum rms drift during that 
time interval of the radio brightness reference point of IM Peg. The rms drift rate was assumed 
to be not larger than one stellar radius over the VLBI observation time span of 8.5 years or 0.06 
mas yr“^ in a and 0.05 mas yr“^ in S. 

The systematic error of the proper motion estimate was found to be mostly comprised of the 
upper limit of the proper motion of Cl in the CRT and of the aforementioned possible maximum 
rms drift of the radio brightness reference point of IM Peg. 

The systematic errors of the parallax and orbit estimates were found to be most clearly defined 
by the position values of B2250+194 relative to those of Cl in 3C 454.3, since the true parallax 
and orbit parameter values should be zero. The same parameter ht used for IM Peg resulted 
in a parallax estimate for B2250+194 relative to Cl of —0.032 ± 0.074 mas, with the statistical 
standard error decreasing to 0.026 mas for = 1- For the five-fold smaller separation of IM Peg 
from 3C 454.3 the systematic uncertainty was estimated to be <0.015 mas in each coordinate. The 
same parameter fit also resulted in orbit parameter estimates for B2250-I-194 relative to Cl much 
smaller than the estimates for IM Peg itself. 

For each of the parameters in Table [3] the total error is also listed. In general it is computed as 
the root-sum-square of the statistical standard error and the estimated systematic error. For the 
position and proper motion values, however, the statistical errors are doubled before computing 
the root-sum-square to allow for correlated noise in the VLBI positions. For the parallax and orbit 
parameter estimates, the systematic errors were too small to contribute appreciably to the total 
standard errors. 


6.3. Comparison of our parameter estimates with previous estimates 


Comparing our results with the most precise previous measurements of pro per inotion and 
parallax, we hnd that our esti mates agree with those in the Hipparcos Catalogue (lESAlll997l ) and 
those of iLestrade et al.l (jl999|) within their respective larger standard errors. Howe ver, we found 
that our es timates slightly disagree with the values in the Hipparcos re-reduction (jvan Leeuwen 
20071 . l2008l i by 1.6 and 2.4 times the combined standard error in and in parallax, respectively. 
This discrepancy in is, however, more than 10 times sma ller than the standard error of either 
the geodetic or the frame-dragging effect (jEveritt et al.ll201lh . and therefore of no consequence for 
the GP-B results. 


Comparing the or bital parameter s , we hnd that our estimate of the time of conjunction is 
consistent with that of lMarsden et al.l (120051 ) and the combination of our estimates of orbital in¬ 
clination and semimajor axis are consistent with their value of asinz. This agreement, in turn, is 
consistent with the reasonable expectation that the radio emission is on average centered on the 
primary of the IM Peg system, and that it orbits with the same inclination. Any offset in phase 
of the radio orbit from that of the primary corresponds to less than one-forth of the radius of the 
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Table 3. Final IM Peg parameter estimates’^ 


Parameter 

Estimate 

Stat. stand, error 

System, error^ 

Total stand, error‘s 

Non-orbit parameters: 

OL at epoch 2005.08^ (errors in mas) 

22** 53“ 2?258612 

0.12 

0.33 

0.40 

5 at epoch 2005.08^ (errors in mas) 

16° 50' 28('16005 

0.13 

0.29 

0.39 

(mas yr~^) 

-20.833 

0.026 

0.073 

0.090 

(mas yr“^) 

-27.267 

0.030 

0.074 

0.095 

Parallax (mas) 

10.370 

0.074 

<0.015 

0.074 

Linear model orbit parameters:^ 

Asa (mas) 

-0.59 

0.10 

<0.1 

0.10 

(mas) 

-0.66 

0.11 

<0.1 

0.11 

Aca (mas) 

0.15 

0.09 

<0.1 

0.09 

Ac 5 (mas) 

-0.23 

0.11 

<0.1 

0.11 

Alternative orbit parameters:® 

Semimajor axis (mas) 

0.89 

0.09 

<0.1 

0.09 

Axial ratio^ 

0.30 

0.13 

<0.1 

0.13 

P.A. of ascending node^ (deg) 

40.5 

8.6 

<8 

8.6 

T^onj (heliocentric 30)-^ 

2450342.56 

0.44 

<0.4 

0.44 


^■Taken from Paper V. 

^See text. 

^The total error computed as the root-sum-square of the statistical standard error and the estimated systematic error, except for the position and 
proper motion parameters where we first doubled the standard errors before computing the root-sum-square. The upper bounds on the systematic 
errors in the orbit terms apply to the mean orbit of the radio emission, and not to the corresponding orbital terms for the stellar binary. For more 
information, see text. 

^The position given is the estimated position of the center of mass of the IM Peg binary at epoch JD 2453403.0 (2005 Feb 1, ~2005.08), the 
approximate midpoint of the GP-B science data. Along with the proper motion, the position is specified in the (J2000.0) coordinate system and 
is closely tied to the ICRF2 iFey et al.||200^b 

^= Ma cos 5. 

^In the linear model, the orbital contribution to IM Peg’s position at time T is Asa sin [27r(T — T’conj)/-^] + Aca cos [27r(T — Tconj)/^] a 
and AsS sin [27r(T - Tconj)/-P] + Ac 5 cos [27r( T — Tcon\)/P\ where P = 24.64877 d is the (fixed) orbital period and Tconj is the (fixed) time 

of conjunction, JD 2450342.905, adopted frorn fMarsden et al.| j2005l) . 

® Computed by iteration to convergence so that the orbit on the sky given by the weighted least-squares estimates of the alternative orbit 
parameters is identical to that given by the weighted least-squares estimates of the linear orbit model parameters. 

^The ratio of the minor axis to the major axis of the sky-projected orbit. 

^See Figure[^for illustration of the orbit geometry. The orbital motion on the sky is counterclockwise. 

•^Time of conjunction for the radio emitting region, that is for the conjunction nearest the one with the primary in back, i.e., at its greatest 
distance from us, for the optical orbit ofiMarsden et al.| l|2005l). 
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primary. 


6.4. Erratic radio emission around the primary of IM Peg 


With position at epoch, proper motion, parallax, orbital motion, and average location relative 
to the center of the primary determined, how can the radio emission be described and where does 
it appear relative to the primary of IM Peg from epoch to epoch? Figure [9] displays a set of VLBI 
images of IM Peg relative to the center of the primary in its orbit about the binary’s barycenter. 


The analysis of the images showed that the radio emission is highly variable and frequently 
partly circularly polarized. The morphology is also variable with the average size of the radio 
emission slightly larger than the disk of the primary. The positions of the peaks of the emission 
regions are scattered over an area on the sky slightly larger than the disk of the primary. The 
scatter has a tendency to be distributed along the orbit normal and expected primary’s spin axis. 
Comparison with simulations suggest that the bri ghtness peaks preferentially occur near the polar 


regio ns similar to the dark spots in the optical (jBerdvugina fc MarsdenI 120061 : iBerdvugina et al 


2000l ). The height of the emission is likely close to the surface of the primary with 2/3 of the 


emission peaks located within 0.25 times the stellar radius above the surface. The radio emission 
may be due to flares linked to a possibly dipolar magnetic field, whose axis is normal to the orbit 
plane, as could be expected giving the tidally locked rotation of the primary. For a movie of the 
star based on these VLBI images, see the website, www.yorku.ca/bartel/impeg.mpg, or the online 
version of Paper VII. 


7. Summary and conclusion 

Our series of 35 VLBI sessions between 1997 and 2005 resulted in the most comprehensive and 
detailed radio investigations ever made of a star. 

1. The proper motion of IM Peg relative to the distant universe is —20.83 ± 0.09 mas yr“^ in 
a and —27.27 ± 0.09 mas yr“^ in 5, where the errors are intended to represent one-standard 
deviation errors and include estimates of systematic errors. 

2. These results met the pre-launch requirements of the GP-B mission of standard errors no 
larger than 0.14 mas yr“^ in each coordinate so as not to discernibly degrade the estimates 
of the geodetic and frame-dragging effects. 

3. The parallax of IM Peg is 10.37 ± 0.07 mas, corresponding to a distance of 96.4 ± 0.7 pc. 

4. The proper-motion and parallax estimates agree with previous estimates, including those in 
the Hipparcos Catalogue, within the combined standard errors, with the exception of the 
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Fig. 9.— A selection of VLBI images of IM Peg at 8.4 GHz with their observing date and FWHM 
of the convolving beam. North is up and east to the left. Both the contours and the color scale 
show the brightness. The contours are drawn at 10, 20, 30, ..., 80, 90, and 98% of the peak 
brightness, starting with the first contour above 3x the rms background; and contours at 50% 
and above are drawn in black. The center of each panel (red dot) is the fit position of the star’s 
center, as derived from the astrometric results in Paper V. In other words, the coordinate origin 
in our radio image should approximately represent the center of the disk of the primary star. Th e 


red c ircle indicates the angular size of the primary star (radius of 13.3 ± 0.6 Rq] iBerdvugina et al 


19991 b The cyan dotted ellipse shows the orbit of the primary as in Paper V and Paper VI. Figure 


taken from Paper VII, where also the complete set of our VLBI images is given. 
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revised Hipparcos results of Ivan Leeuwenl (|200a) where the agreement is within 1.6 and 2.4 
times the respective combined standard errors. Any discrepancy at these levels is however of 
no conseqnence for the GP-B results. 


5. The radio emission of IM Peg is highly variable. 

6. The VLBI images of IM Peg show emission regions slightly larger than the disk of the pri¬ 
mary and their peaks scattered about a similarly sized area, consistent with being located 
preferentially near the polar regions and within ~0.25 stellar radii above the snrface. 

7. The VLBI images of IM Peg were assembled for a movie of the star (see the website of the 
hrst author, www.yorku.ca/bartel/impeg.mpg, or the online version of Paper VII). 
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